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Abstract. We propose an experimental implementation of a topological superfluid 
with ultracold fermionic atoms. An optical superlattice is used to juxtapose a ID gas 
of fermionic atoms and a 2D conventional superfluid of condensed Feshbach molecules. 
The latter acts as a Cooper pair reservoir and effectively induces a superfluid gap in 
the ID system. Combined with a spin-dependent optical lattice along the ID tube and 
laser-induced atom tunneling, we obtain a topological superfluid phase. In the regime 
of weak couplings to the molecular field and for a uniform gas the atomic system is 
equivalent to Kitaev's model of a p-wave superfluid. Using a numerical calculation we 
show that the topological superfluidity is robust beyond the perturbative limit and in 
the presence of a harmonic trap. Finally we describe how to investigate some physical 
properties of the Majorana fermions located at the topological superfluid boundaries. 
In particular we discuss how to prepare and detect a given Majorana edge state. 
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1. Introduction 

Since the discovery of topological insulators [TJ [2] the search for exotic states of matter 
with non-trivial topology such as a topological superconductor has generated strong 
interest. A non-trivial topological phase is characterized by a gap in the bulk elementary 
excitations and a non-local topological order associated with gapless states at the system 
edges, defects or vortex cores [3] E] EJ El [JJ [8]. In topological superconductors, these 
edge states can be described as Majorana fermions, which are quasi-particles equal to 
their own antiparticles [9]. Besides the fundamental interest in observing Majorana 
fermions a great motivation for their study relies on their non-Abelian statistics, which 
could be used for topological quantum computation immune to decoherence from local 
perturbations [TO] QH Q21 Q21 Q3] . 

The simplest model of a topological superfluid was introduced by Kitaev in 2001 
[6]. It consists in a ID lattice system of spinless fermions with a p-wave superfluid gap, 
described by the Hamiltonian 

#Kitaev = J2 ~\)~ J ( £ ^+l + h - C + A (^+1 + hX -) ' (1) 

i ' 

where cj creates a fermion into the lattice site i, and /x, J, A respectively represent 
the chemical potential, the tunneling amplitude and the superfluid gap. However 
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this 'toy model' cannot be directly realized with solid-state systems. Subsequently 
several proposals described how to create topological superconductors in solid-state 
heterostructures combining superconducting proximity effect with a topological insulator 
or a semiconductor with strong spin-orbit coupling [151 EE EH HI HH]- These proposals 
led to the first observations in 2012 of boundary quasiparticles behaving as Majorana 
fermions [201 ED E2 E3] . 

Ultracold atomic gases constitute an alternative system for realizing a topological 
superfluid [211 EHl ES EH EH EH] • The extreme purity of these gases could allow one to 
realize model systems in a clean and controlled manner [30 J . As described in [3lT 132] 
the study of atomic Majorana fermions would also benefit from measurement techniques 
specific to atomic physics. While inducing topological superfluidity using resonant p- 
wave interactions cannot be achieved due to large inelastic loss rates [211 E31 El], one 
expects that, in a quasi-2D Fermi gas, combining s-wave interactions with Rashba-type 
spin-orbit coupling may lead to the formation of a p x + ip y superfluid, which has a 
non-trivial topology [35l EH [371 EHl EHJ HO] . Alternatively, by analogy with proximity- 
induced superconductivity [41j, superfluidity could be induced by the coherent coupling 
to a molecular Bose-Einstein condensate rather than by interatomic interactions [28] . 

In this article we propose an experimental implementation of a topological 
superfluid with ultracold Fermi gases based on proximity-induced superfluidity [28]. We 
use a superlattice potential to juxtapose a quasi-2D gas of condensed Feshbach molecules 
and a set of ID tubes. The coupling of atom pairs between both systems effectively 
creates an s-wave superfluid gap in the ID tubes. In addition, the atoms in the ID 
tubes are held in a spin-dependent optical lattice, in which atom tunneling is induced 
by the appropriate laser transitions. The scheme makes use of several ingredients - such 
as quasi-2D Fermi gases [121 HHJ SU 05] , spin-dependent potentials [36l Wf[ I3K1 1391 150] 
or laser-induced atom tunneling [511 [52] - that were already developed in the cold atom 
community. We show that this setup leads to the realization of a topological superfluid, 
which is equivalent to Kitaev's model in the perturbative regime of weak coupling to the 
molecular field. For stronger couplings we obtain a robust topological superfluid with a 
large bulk gap value, both for homogeneous and harmonically trapped systems. Finally 
we describe some measurable physical properties of a trapped atomic gas in this lattice 
geometry, which should exhibit a topological superfluid phase with Majorana fermions 
at its boundaries. In particular we discuss how to prepare and detect a given Majorana 
edge state. 

2. Proximity-induced superfluidity 

Before describing the implementation of Kitaev's model, we first discuss the realization 
of proximity-induced s-wave superfluidity in a two-dimensional geometry. We consider 
a pair of quasi-2D atomic systems (planes A and B) that can be created using a double- 
well potential along z (see figure [T]d). The plane B is filled with a two-component Fermi 
gas, with effective spins labelled |t) and We assume that a Feshbach resonance is 
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used to create weakly-bound Feshbach molecules [HJ EH], and that the gas is cooled 
sufficiently to form a nearly pure molecular Bose-Einstein condensate (mBEC). When 
coupled to the plane A as described below, this system will act as a reservoir of Cooper 
pairs and effectively induce a superffuid gap in A [28J. 



2.1. Effective Hamiltonian of a gas coupled to a molecular BEC 

The reservoir of condensed molecules in plane B is described by a coherent state with a 
mean number of molecules N. The 2D wavefunction of a single molecule is written as 

</>mol(Pt> Pi) = ^0rclOf - Pi), 

where L is the size of the quantification box along x and y and 

W(p) = —^K (p/a) 
'Tea 



is the 2D molecular wavefunction for the relative motion. Here a denotes the 2D 
scattering length and Kq is a modified Bessel function of the second kind. The molecular 
binding energy will be written as Ej, = H 2 /ma 2 . 

The coupling from the reservoir to the plane A occurs via atom tunneling through 
the A — B barrier, described by a Hamiltonian 

H z = -J z J dp (V^ )t (p)^ B)t (p) +$A, i {p)^ B ,i{p)) + h - c ; 

where J z is the single-atom tunneling amplitude and ffi x CT (p) denotes the creation 
operator of a particle of spin a in the plane X. 

We expect that tunneling of atoms towards the plane A leads to the dissociation of 
molecules, leading to pairs of atoms of opposite spin in both planes A and £(U These 
dissociated atom pairs can be described by an effective Hamiltonian using a semi-classical 
approximation, which consists in replacing the mBEC by a classical field. By analogy 
with the Glauber unitary transform for coherent states of the electromagnetic field [53] 
we define the effective Hamiltonian as 

Hz,eft — -^mBEC ^ mBEC ' 

where 

AnBEC = e- v ^ /2 exp(v / iV(^ i - VW)) 
is the displacement operator that brings the mBEC to the vacuum state and 

^mol = J d Pi d Pi 0mol(Pt, P4^B,tM$uW 

creates a molecule in plane B. 

The effective Hamiltonian can be calculated explicitely as a series in kpa, where 



kjp — 2y / 7rviV/L is the Fermi momentum in the plane B (see Appendix A). As we are 

| Note that the tunneling of bound molecules to the plane A is prevented by the spin-dependent optical 
lattice described in the following. 
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Figure 1. a. General scheme of the setup, consisting in quasi-lD gases of fermionic 
atoms coupled to a 2D gas of condensed molecules. The tunneling of atoms pairs into 
the ID tubes acts as a superfiuid gap applied externally. A spin-dependent optical 
lattice is applied along the tube, leading to an effective lattice with spin f (I) particles 
localized in the even (odd) sites. Atom tunneling inside the tube is driven by Raman 
lasers, leading to the realization of Kitaev's model of a p-wave superfiuid. b. Scheme 
of the atom pair tunneling from A to B, resonant for = 2 5. 



interested in the molecular regime kpa <C 1 we only keep for simplicity the first term in 
/cpa, leading to a Hamiltonian H ZtB s = H z + H p &hs, where the additional term 

^pairs = - J *7^ J d Pt d Pl 0rel(Pt ~ Pi) (^B.tW^U W + ^(P^B^Pi) + h - C - 

describes the creation (or annihilation) of pairs of particles dissociated from the mBEC 
We show in the |Appendix A| that the neglected terms of the series have minor effects 
for the parameters used in the following. 

An additional energy offset 5 is introduced in the plane B in order to match the 
resonance condition 25 = E^ (see figure [l]b), described by the Hamiltonian 

H 5 = 5 J dp (^, t (p)^ At (p) + -0l )i (p)^(p)) ■ 
The total effective Hamiltonian is then given by the expression 

Ha = H z + impairs + Hg. (2) 



2.2. Superfiuid gap in the weak coupling regime 

A simple picture of the system can be obtained in the limit J z ^ 5 where at low energy 
all non-condensed atoms are located in the plane A and the plane B is only virtually 
populated by single atoms. In second-order perturbation theory, the system can be 
described by a Hamiltonian restricted to the plane A 

rr(2) _ p (H z + impairs) 2 p 
n A — "A ^ "At 

where Pa projects on the Hilbert space of wavefunctions localized in the plane A. A 
straightforward calculation leads to 

#i 2) = [dp t dp i A( Pt -p i )^(p t )^ Al (p i ) + h. C ., (3) 
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A(p)= - ^-ffa(p), 
V 71 " <-> 

which describes a system of two-component fermions with an order parameter A(p). 
The latter has an s-wave symmetry inherited from the symmetry of the molecular 
wavefunction. 



3. Creating a topological superfluid 

The technique of proximity-induced superfluid gap presented in section [2] can be used to 
create a topological superfluid. For simplicity in this section we restrict the discussion 
to the limit of perturbative coupling to the mBEC for which the atom dynamics is 
restricted to the plane A. 

Optical lattices along x and y are applied in the plane A only. The y lattice 
freezes the motion along y, leading to a quasi-lD geometry. The x-lattice is a spin- 
dependent tight-binding optical lattice as described in figure [T^,: each spin state is held 
in a sublattice of period 2d, with the two sublattices shifted by a distance d. This system 
can be viewed as a lattice of spacing d with spin f (J,) particles localized in the even 
(odd) sites 2% (2i + 1, respectively). The practical realization of this lattice structure is 
described in section 15.21 

The physical system restricted to the lowest band of the optical lattice can formally 
be viewed as a spin-polarized system of fermions on a discrete lattice. Expanding the 
fermionic field operators on the basis of Wannier functions: 



^m(p) = ^2 w x(x - (2i + l)d)w y (y)c. 



t 

2i+l' 



we rewrite the Hamiltonian (|3j) as 

H {2) - V A (2) r f r f + h c 



where 



A 2,2i+1 = J d P^ d Pl ^rcl(Pt - Pi) >< 

Wy(y t )w y (yi)w x (x t - 2id)w x (x i - (2j + l)d). 

i (2~) 

An analytic expression for A 4 J can be obtained assuming that the 1/e size a x of the 
Wannier function along x is much smaller than the lattice spacing d, and that the extent 
of the wavefunction along y is much larger than the 2D scattering length a: 
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In the following we will use values for the scattering length a = d/2, for which 
the Aj j is strongly suppressed except for nearest-neighbor pairs of sites, leading to the 
Hamiltonian 

= J2 4(4+1 + 4-i) + h-c = A( 2 ) + h-c, 

j i 

where 

A (2) = Ag) +1 = -2^ J -fk F a x e- d l\ (4) 

This Hamiltonian differs from the gap term in Kitaev's model due to the (—1)* factors. 
However, the latter can be eliminated by redefining the operators q according to 
c 2 i -> (-l) l d 2 i and c 2i+ i -> {-l) % d 2 i+i- 
The complete Kitaev model 

-^Kitacv = ^2 ~^ y^di ~ 2) ~ J (d\d i+1 + h.c.j + A (d\d\ +l + h.c.) 

i 

is then obtained by introducing a small additional energy offset \x of the plane A 
and inducing the required atom tunneling J between neighboring sites. The latter 
is performed using Raman transitions that flip the atom spin, which induces atom 
tunneling with the required phase (—1)* of the hopping operator d\d i+ i = (— l)*cjc i+ i 
expressed with the physical operators C; [51J. The required laser setup is described in 
section 15.21 



4. Bogoliubov spectra and Majorana fermions 

The topological character of a topological superfluid can be characterized by its spectrum 
of Bogoliubov excitations. Let us remind the main properties of Bogoliubov excitations 
in the case of Kitaev's model, described by the Hamiltonian (|T| [6]. In the case of a 
homogeneous ID chain with periodic boundary conditions, the spectrum of Bogoliubov 
excitations is gapped when A > 0, except when the chemical potential is fine-tuned to 
\i = —2 J or 2 J. These points indicate phase transitions between topologically different 
systems. The phase fi < —2J can be smoothly connected to the trivial vacuum by 
decreasing \i to —00 without closing the gap, it is thus non-topological. The phase 
H > 2 J is essentially equivalent to /x < —2.1 since both phases are related by a particle- 
hole transform. On the other hand, the phase —2 J < /j < 2J has a non-trivial topology 
which can be characterized by a Z 2 topological invariant [6]. 

Similarly to edge states in the topological quantum Hall phases, for a chain with 
open boundary conditions the non-trivial topology is associated with boundary states. In 
the case of ID topological superfluids the boundary states form a pair of Majorana modes 
occuring at zero energy, i.e. in the middle of the gap of bulk Bogoliubov excitations. 
These states are highly non-local - delocalized at both ends of the topological superfluid 
- and become exponentially degenerate when increasing the size of the topological phase. 

We present in this section spectra of Bogoliubov excitations calculated for realistic 
experimental parameters, using the Hamiltonian described in sections [2] and [3] (without 
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Figure 2. Amplitude of the bulk superfluid gap Abuik as a function of the tunnel 
coupling J and the chemical potential [i calculated for a system of 30 lattice sites with 
periodic boundary conditions. The topological transition corresponds to the lines for 
which the gap cancels; their location agrees well with the prediction of Kitaev's model 
/i = ±2 J (dashed white lines). 



performing the perturbation theory to second order in J z /8). The system parameters 
are the following: the lattice spacing d is equal to half the wavelength A of the optical 
lattice (see section 5.2), and we choose a 2D scattering length value a = d/2 so that 
the proximity-induced gap is restricted to nearest neighbours. This scattering length 
value corresponds to a molecular binding energy E b ~ 0.8 E r , where E r = h 2 /2m\ 2 
denotes the recoil energy. The molecule density in B is fixed to a value n mo i = 0.16 d~ 2 
corresponding to a Fermi energy Ep = E^/A, i.e. on the BEC side of the BEC-BCS 
crossover of s-wave Fermi superfluids. This condition ensures that at low temperature 
the condensed fraction of the mBEC tends to unity [51] . We choose an x-lattice depth 
value V x = 5E r , which leads to a spontaneous atom tunneling of negligible amplitude 
6.2 x 10~ 4 E r . The band gap is then equal to 2.0 E r , a value larger than the energy offset 
5 = E b /2 = 0AE r , as required by the single-band assumption made above. 

For the Raman-induced tunneling rate we take a value on the order of the p-wave 
gap A^ 2 ) calculated in second-order perturbation theory, which should maximize the 
system bulk gap Abuik in the perturbative limit J z <C 5 (more precisely Abuik = 2 A for 
J = A for Kitaev's model [6]). 



4-1. Homogeneous topological superfluid 

We discuss in this section the case of a homogeneous system. 

We first calculate the spectrum of Bogoliubov excitations for a homogeneous system 
of 30 lattice sites with periodic boundary conditions. It exhibits a bulk gap whose 
amplitude Abuik is plotted in figure [2] as a function of the tunnel coupling J and the 
chemical potential \i. We observe a cancellation of the bulk gap on two lines of the 
phase diagram ji ~ ±2 J, in agreement with the prediction of Kitaev's model. The lines 
indicate the location of the topological transition. 
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Figure 3. a. Spectrum of Bogoliubov excitations (red dots) calculated for J z = 0.2 5 
and J = A^ 2 \ and compared with the prediction of Kitaev's model with J = A = A' 2 ) 
(black dots), b. Density distribution along x of a zero-energy Majorana state, in the 
planes A (red line) and B (blue line), revealing the non-local character of Majorana 
states. In the perturbative regime J <C 5 the population in B remains small. 



We then repeat the calculation of Bogoliubov spectra for homogeneous systems 
with sharp boundaries, for which zero-energy Majorana end states are expected in the 
topological phase. We first consider a small tunneling rate value J z = 0.2 5. As shown in 
figure [3^, the calculated spectrum of Bogoliubov excitations is very similar to the one of 
Kitaev's lattice model, with a bulk gap and two Majorana states at zero energy. The full 
excitation spectrum quantitatively differs from Kitaev's model due to additional effects 
such as the creation of atom pairs in non-neighbouring lattice sites; yet the bulk gap 
value is in good agreement with the second-order perturbation result of equation Q. 

When increasing the tunneling rate value beyond the weak coupling regime J z <C 5, 
we find that the structure of the Bogoliubov spectra qualitatively remains the same for 
tunneling rates J z ~ 5 (see figure |4^i). This should allow one to create p-wave superfiuids 
with large bulk gap values A bulk ~ 0.2 E r , which corresponds to a gap of k B x 40 nK for 
the physical system described in section [5] (see figure As the gas temperature has to 
be much smaller than the Fermi energy Ep ~ 0.2 E r in order to create a molecular BEC 
in plane B [55] , we expect the p-wave superfluid to be essentially in the zero-temperature 
regime for J z ~ 5. 

The density probability of the Majorana states at zero energy is shown in figure [lb 
in the case J z = 5. In the plane A essentially only the lattice sites at the edge of the 
system are populated. We note that in that case the probability to occupy the plane 
B is not negligible (about 30%). The total atom density probability is almost uniform 
and close to half filling (see figure |4]i). 

4-2. Topological superfluid in a harmonic trap 

Homogeneous systems with sharp boundaries are difficult to realize with ultracold atom 
systems. We thus now focus on gases held in a harmonic trap V(x) = ^mu 2 x 2 . In a 
local density approximation picture, the local properties of the gas are expected to be 
close to the ones of a homogeneous gas with a local chemical potential /i — V(x). A 
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Figure 4. a. Spectrum of Bogoliubov excitations (red dots) for a homogeneous 
system with sharp boundaries, calculated for J z — 5, and J = 2A^ 2 K It exhibits a 
bulk gap Abuik = 0.18 E T and a pair of zero-energy Majorana states with a residual 
splitting A s ~ 10~ 12 _E r . b. Evolution of the bulk gap amplitude Abuik as a function 
of the ratio J x /6 (red line), for J = A (2) , with A^ 2 ) given by eq. Q. The black 
line represents the prediction of second-order perturbation theory Abuik = 2 A^ 2 '. c. 
Density distribution along a; of a zero-energy Majorana state, in the planes A (red 
line) and B (blue line, offset for clarity). In the strong coupling regime J ~ 5 the 
population in B is not negligible, d. Total density distribution along x calculated 
at zero temperature. Majorana states are not visible in this almost uniform density 
profile. 



gas of chemical potential — 2 J < // < 2 J should thus exhibit a phase separation with a 
topological superfluid phase in the trap bottom (/j l — V(x) > —2 J) and trivial superfluids 
at the edges (p — V(x) < —2 J). We also expect the presence of Majorana zero-energy 
states at the boundary of the topological superfluid. 

For the parameters given in the caption of figure [5j the calculated spectrum of 
Bogoliubov excitations is gapped with a pair of quasi-degenerate Majorana end states 
at zero energy (see figure [5^,). These states are non-local, being the superposition of 
two localized wavefunctions (1/e size of ~ 2.5 d) separated by ~ 21 d (see figure ^jp)- 
The position of these zero-energy modes agrees well with the position of the topological 
transition given by the local density approximation. Note that the Majorana states are 
not directly visible in the total density profile; yet at the location of the topological 
transition the gas compressibility sharply varies (see figure [5]i), which could be useful 
to establish the phase diagram of the system [561 EZ] • 

A figure of merit of the topological protection can be given by the residual splitting 
between the two Majorana states. As shown in figure [5]d, the topological protection 
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Figure 5. a. Spectrum of Bogoliubov excitations (red dots) for a harmonically 
trapped system, calculated for J z — 1.5 6, J = 2 A^ 2 \ u) = 2n x 200 Hz. It exhibits a 
bulk gap Abuik = 0.22 E r and a pair of zero-energy Majorana states with a residual 
splitting A s ~ lO -10 ^. b. Evolution of the energy splitting A s between Majorana 
states (red dots) and of the bulk gap Abuik (orange dots) as a function fo the trapping 
frequency u>. While the bulk gap essentially does not depend on w, the splitting between 
Majorana states strongly decreases when decreasing the w value, i.e. increasing the 
separation between Majorana fermions. c. Density distribution along a; of a zero- 
energy Majorana state, in the planes A (red line) and B (blue line), d. Total density 
distribution along x calculated at zero temperature. 

improves by increasing the separation between the Majorana fermions, which can be 
performed by opening the trap along x. 

4-3. Probing the Majorana fermions 

As shown in figure [5^,d, the presence of Majorana end states cannot be revealed from 
the total atom density. Photo-emission spectroscopy [58], which directly addresses the 
structure of the Bogoliubov excitation spectrum, is a convenient probe for these states 
[3TL [32J. It consists here in coherently coupling the spin ||) atoms to an unpopulated 
internal state \a), which can be probed either in situ or after time-of-flight. The 
photoemission performed at a given frequency Q gives rise to a transfer of atoms 
occupying the Bogoliubov excitations of energy E matching the resonance condition 
HQ = hu>o — E, where Uq is the resonance frequency for a single atom. The energy of 
the final state |a) is neglected for simplicity. 
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Figure 6. Position-resolved photoemission spectrum calculated with the parameters 
of figure [5j The Majorana states probed at zero energy are localized at the edges of 
the topological superfluid. 



4-3.1. Density probability of Majorana states. In order to probe the density probability 
of Majorana states, we probe the target state density distribution right after the 
photoemission pulse. The target spin state \a) can be chosen as another Zeeman 
state that would be trapped in x lattice, with a photoemission coupling provided by 
a radiofrequency field. We assume that the photoemission spectroscopy is performed 
in the linear perturbative regime, and that the target spin state |a) can be imaged in 
situ with single-atom sensitivity [59, 60J. The density probability of \a) atoms as a 
function of the driving frequency is shown in figure |6j The bulk Bogoliubov excitations 
are associated with a broad energy spectrum above the superfluid gap. The Majorana 
boundary states correspond to a narrow feature at zero energy, with photoemitted atoms 
located at both ends of the topological superfluid phase. 

4-3.2. Momentum distribution of Majorana states. The photoemission technique can 
also be applied to extract the momentum distribution of Majorana states, by imaging 
the target state \a) after time-of-flight. 

For this measurement the target state |a) should be untrapped so that its 
eigenstates have a well-defined momentum. A possible choice for |a) is a long-lived 
excited state^J for which the optical dipole force of wavelength A used for the optical 
lattice would be very weak [62]. The photoemission coupling would be achieved using 
the ultranarrow optical transition between the ground electronic state and the state 

The expected photoemission spectrum (shown in figure [7^,) can be divided in a 
broad spectrum of bulk excitations and a narrow feature at zero energy corresponding 
to the Majorana states. The oscillatory behavior at zero energy reflects the quantum 
superposition of the wavefunction of Majorana states at both ends of the topological 
superfluid. Besides demonstrating the non-local character of the Majorana states, this 
probe also allows one to discriminate between the two zero-energy states, whose phase 
of the interference pattern differ (shown in figure [Tjo). 

§ such as the 4f 9 ( 6 H°)5d6s 2 7 H^ metastable state of Dy (1/e lifetime of 7 ms), coupled to the ground 
electronic state by an optical transition at 1322 nm [61] . 
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Figure 7. a. Momentum-resolved photoemission spectrum, exhibiting an oscillatory 
behavior at zero energy reflecting the delocalization of Majorana states at both ends of 
the topological superfluid. The small rectangle at zero energy represents the spectrum 
of Majorana states shown in b. b. Momentum-resolved photoemission spectrum rii(k) 
for both Majorana states i = 1,2, differing from the phase of the fast oscillations. 
The slow modulation (period ir/d) reflects the wavefunction modulation by the lattice 
potential, while the fast oscillation is due to the delocalization of the wavefunction at 
both ends of the topological superfluid. 



4-3.3. Adiabatic preparation of a given Majorana state. In the previous section we 
described how to discriminate between Majorana states using momentum-resolved 
photoemission spectroscopy. This raises the question of the preparation of the system in 
a given Majorana state |31j . which is an important prerequisite to the manipulation of 
Majorana fermions for quantum information purposes. We propose to prepare a given 
Majorana state by crossing adiabatically a phase transition from a trivial to a topological 
superfluid phase. This can be achieved by changing adiabatically the chemical potential 
fi, which can be performed by dynamically varying the double- well energy offset (see 
section 5.2.1 ). 

In the local density approximation framework, one expects a topological superfluid 
phase to nucleate in the trap bottom when the chemical potential corresponds to the 
critical value for topological superfluidity /i c at the bottom of the trap (fi c ~ —2 J). As 
shown in figure |8^, the nucleation of the topological superfluid is associated with the 
creation of Majorana zero-energy states from two bulk states (one from the hole sector, 
one from the particle sector) when increasing the chemical potential value across /i = fi c . 
By increasing adiabatically the chemical potential value from a trivial to a topological 
phase, one thus expects the system to populate the Majorana state which is adiabatically 
connected to the hole sector of the trivial superfluid. In order to test the feasability of 
this scheme, we simulated the effect of the finite duration of the chemical potential ramp 
by solving the Schrodinger equation for our system. As shown in figure |HJd, fidelities 
J 7 > 0.9 can be achieved with the parameters described in the caption of figure [8] for 
ramp durations of about 120 ms. Note that shorter adiabatic paths could be used if 
one crosses the topological transition with a lower atom number in plane A (i.e. higher 
trapping frequency u), and then increase the atom number in A by opening the trap 
far from the topological transition so that Majorana states are deeply topologically 
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Figure 8. a. Energy of Bogoliubov excitations as a function of the chemical 
potential [i. The parameters are the same than for figure [5] expect with a trap 
frequency uj/2t: — 600 Hz. The bulk hole sector is pictured in blue. Starting with 
a system prepared at low temperature in the trivial phase (ji ~ —2E r ) the Majorana 
state connected to the hole sector of the trivial superfluid (in red) can be prepared 
by increasing adiabatically the chemical potential up to /i = 0. b. Expected 
fidelity J 7 as a function of the ramp duration T (black dots). The solid line is an 
exponential fit to the numerical data. A fidelity J- > 0.9 is expected for ramp durations 
T > 3300 K/E r ~ 120 ms. 



protected. 

5. Description of the experimental setup 

5.1. Choice of the atomic species 

The proposal described above makes uses of a spin-dependent optical lattice potential. 
For fermionic alkali atoms ( 6 Li and 40 K) spin-dependent optical dipole potentials are 
associated with large heating rates due to spontaneous emission [361 EH], Decause of 
the small fine-structure splitting of the P electronic levels compared to the transition 
linewidth. In order to circumvent this issue we propose to use a fermionic atom of 
the Lanthanide family, such as 167 Er or 161 Dy [631 EU [65] . These atoms benefit from 
narrow optical transitions that can be used to create spin-dependent dipole potentials 
with small heating rates. For example, strong spin-dependent potentials can be created 
for Dy atoms using optical lattices at A = 530 nm, close to the 530.3 nm transition 
(electronic angular momenta J — 8 — > J' = 7, linewidth r = 2tt x 184 kHz). 

Due to the large magnetic moment of these atoms, we expect large dipolar relaxation 
rates, except for mixtures of the two Zeeman states of lowest energy in the case of a 
fermionic species |66j . As shown in the case of 168 Er Bose-Einstein condensates, the very 
large number of Zeeman states in the electronic ground state should lead to numerous 
Feshbach resonances, making it possible to create a molecular Bose-Einstein condensate 
for large scattering length values [65J. 
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Figure 9. a. Scheme of the Dichromatic superlattice potential along z (black line), 
made of the superposition of two standing waves of wavelength A (blue line) and 2 A 
(red line). The energy offset 8 in the resulting double- well potential can be controlled 
using a relative offset in position of the two standing waves, b. Amplitude of the 
standing wave along z used for generating the y lattice, which vanishes at Zg. c. 
Superlattice potential with the Wannier function in plane B centered on zg. 

5.2. Optical lattice configuration 

We describe in this section the laser configuration required to create the double-well 
potential along z, the plane-dependent lattices along x and y, and the Raman lasers 
driving atom tunneling in the x lattice. 

5.2.1. Superlattice potential along z. The atoms are confined in a bilayer system of 
two planes A and B using a double-well potential. The latter can be created using 
a bichromatic optical lattice along z, i.e. the superposition of two standing waves of 
wavelength A and 2A, leading to a potential 

V S l(z) = V 1 cos 2 \kz) + V 2 cos 2 (kz/2 + 0), (5) 

where k = 2ix/\ [67] (see figure (9^,). Atom tunneling towards other planes is suppressed 
using lattice depths of about 20 E r . The lattice intensities V\, V 2 and relative phase </> 
enable one to control independently the tunnel coupling J z between A and B and the 
value of the energy offsets 8 and \i. 

5.2.2. Plane- dependent lattices along x and y. An optical lattice is used to freeze the 
atom motion along y in plane A. The atom motion in plane B, which contains the two- 
dimensional niBEC, should remain unaffected so that phase coherence in the mBEC is 
maintained. 

The plane-dependent optical lattice along y can be created using the interference of 
a travelling wave propagating along u = cos a e y + sin a e z with a standing wave along 




A/2 



Figure 10. a. Scheme of the electronic states used for creating spin-dependent 
optical lattices with a J- > J' = J — 1 optical transition, b. Lattice potential along x 
obtained by interfering a standing wave along z with a plane wave along x. The required 
spin-dependent lattice configuration is achieved for specific choices of polarization, c. 
Geometry of the complete laser configuration generating the optical lattice potential 
and the laser-induced tunneling in the x-lattice. 

z (see figure [lOfc), creating a potential 

V (r) oc | E u e iku *- r + E z cos{kz + <j)') \ 2 

= \E Z \ 2 cos 2 {kz + 4>') + 

2\E Z E U \ cos(kz + (J)') cos (k( cos ay + sin a z)) + \E U \ 2 . 

The first term represents a simple optical lattice along z that contributes to the 
superlattice potential. The second term consists in an optical lattice along y of lattice 
spacing d y — A/ cos a whose amplitude depends on the plane positions za and z B . As the 
z standing wave has a node on the retroreflection mirror, the phase <j)' can be controlled 
with high precision by varying the laser frequency over typically a GHz range J67J 68J . 
This allows one to choose the value of 0' so that the y lattice amplitude vanishes at zb, 
leaving the motion of molecules in the plane B unaffected (see figure [9]d,c). 

Choosing the value a = n/3 corresponds to a y lattice spacing d y = 2 A. For such 
a value d y ^> a the amplitude for atom pair tunneling into different tubes is about 1CT 3 
smaller than for pair tunneling into the same tube. This allows one to consider that the 
tubes are independent from each other. 

The same technique can be used to create the lattice potential along x. As the 
optical lattice wavelength is chosen close to a narrow transition, vector light shifts can 
be used to create the spin-dependent potential in the desired configuration shown in 
figure [10)3 . We checked that using specific laser polarizations we can create the required 
spin-independent lattices along y and z and the spin-dependent lattice along x (see 
Appendix B). 



5.2.3. Laser-induced tunneling In order to realize Kitaev's model atom tunneling in 
the x lattice should be described by the Hamiltonian — l)*c|c i+ i. The required 
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alternation of sign of the tunnel amplitude can be achieved by inducing atom tunneling 
via Raman transitions in the following geometry (see figure [Tofc): one of the Raman 
beams would be the z lattice beam used for generating the x and y lattices, and the 
second beam would propagate along x, with a frequency detuning equal to the |f) -H- \i) 
transition frequency. The hopping matrix element between the sites 2i and 2i + 1 is 
then given by 

hn J dr cos(fc2 + <p z )e ikx W 2i (r)*W 2i±1 (r), 

where W^(r) = w x (x — id)w y (y)w z (z) is the Wannier fonction of the site i of the x lattice 
and Qq is the Rabi frequency of the Raman transition in vacuum. An explicit calculation 
leads to the result 

J2i,2i±i = ±J, J = J dz cos(kz + 4> z )\w z (z)\ 2 J dxe lkx w x (x)*w x (x - 
which is the desired configuration, assuming the gauge is chosen so that J is real. 



2i,2i±l 



6. Conclusions and perspectives 

In this article we have discussed how to realize ID topological superfluids with ultracold 
Fermi gases. We described how to probe the Majorana fermions located at the edges 
of the topological superfluid phase. Further developments could address the robustness 
of the Majorana degeneracy with respect to disorder (691 EO] and atomic interactions 
between neighbouring sites [7H [721 E2 HI] i which could be induced by the large dipole- 
dipole interactions for Dy or Er atom^jj] On a longer term the manipulation of Majorana 
fermions would allow one to exhibit their non-Abelian statistics, which is the basic 
ingredient for developing topological quantum computation [101 HH EG] • The realization 
of braiding operations requires either implementing a network of coupled ID tubes [75] or 
creating a two-dimensional p x +ip y topological superfluid with Majorana fermions bound 
to half-quantum vortices (U [121 E>]- The experimental scheme presented in this article 
could be adapted to a 2D geometry to realize a lattice p x + ip y topological superfluid. 
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Appendix A. Semi-classical approximation 

In this appendix we provide details on the calculation of the proximity-induced gap 
using a semi-classical approximation. For replacing the mBEC coherent state by 

|| For two Dy atoms separated by A/2 = 265 nm the dipole-dipole interaction varies between — 0.03EV 
and 0.016 E r , depending on the magnetic field direction. 
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a classical field we use a unitary displacement of the Hamiltonian by the operator 
AnBEC = e _v ^ /2 exp(A/iV(^ ol - Vmoi))- For the calculation of D^^H^A 
use the formula 



Vbec we 



1 
2! 



A. 



A,B 



e A Be~ A = B + A,B 
This leads to an effective Hamiltonian 



+ 



3! 



A. 



A. 



A, B 



+ . 



J z J dpdp' f(p' - p') ($i jt (p)4,iV) + ^m(p)^b,i(p') + h.c. 
- J z J dp t dp 4 #(p t - p t ) (^B,-\{P\)$A,\.{Pi) + A,t(Pt)^kl(^) + h - ( 



where the first terms describes atom tunneling between the planes A and B and the 
second terms corresponds to the creation/annihilation of atom pairs from the mBEC. 
The functions f(p) and g(p) are given by 



M = E 



l) n+1 /fc F a 



ra=0 



(2n + 1)! 



2tt / 

2n+l 



2tt 



T < r *2n+1 



where i^Q* refers to the function K convolved with itself % times. 

In the main text we only kept for both functions the first term of the series. We 
checked that for the value Ep = E^/A used for the calculations this assumption is 
justified (error on g(p) smaller than 10%). 



Appendix B. Polarization of the optical lattices 



We provide in table |B1| a possible choice for the polarization of the lasers used for the 
spin-independent lattices along z and y, and the spin-dependent lattice along x. 
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